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O ■ Abstract 

■^ ■ 

G> ' The response functions for the extended second and third random phase approximation are 

compared. A second order perturbation calculation shows that the first-order amplitude for 
, ^ . the direct 3p3/i excitation from the ground state cancels with those that are engendered by the 
3 ■ lpl/i-3p3/i coupling. As a consequence nonvanishing 3p3/i effects to the Iplh response involve 



off energy shell renormalization only. On shell 3p3/i processes are absent. 



tPellow of the CONICET from Argentina. 



Many efforts have been devoted during the last few years in developing generalized random 
phase approximations (RPA), which go beyond the standard one-particle - one-hole {Iplh) 
approach [|I|. This has been accomplished by including additional correlation effects in both 
the ground state and the excited states § |, |, |, |, 0, |, |10|, |lll |12|, |13|, |14|, |15|, [TBI. The 
reasons for that were mainly: i) the problem of the missing strength in the Gamow- Teller 



(GT) resonances, induced by {p,n) reactions ^^ |l^, and ii) the issue of the missing charge 
and missing dip-strength in quasielastic electron scattering [|l^. In particular, the extended 
second RPA (ESRPA), which explicitly includes the 2p2h ground state correlations (GSC), was 
extensively used to describe the above mentioned nuclear excitations 0, |, |, |, [l^, |Tl|, [12|, |13|, 
0. Yet, it is self evident that when the 2p2h admixtures are present in the ground state, the 
external excitation field can lead, not only to the Iplh and 2p2h states in the final nucleus, but 
also to the 3p3h states. However, as the ESRPA does not involve the 3p3h propagator these 
excitations cannot appear within the response function as real on the energy-shell processes. 
Recently the 3p3h degrees of freedom were explicitly included within a Tamm-Dancoff approach 
(TDA), and their effects on the non-energy- weighted GT sum-rule were discussed |T^. Also an 
extended third RPA (ETRPA), which possesses as the TDA limit the formalism developed in 
ref. [^, has been used to study the effects of 3p3h excitations on the static strength function 
for quasielastic electron scattering |TB . 



The purpose of this paper is to present some results for on the energy-shell 3p3h effects 
in the response function. This is done in the context of the full ETRPA approach which is 
therefore reviewed below. The nature of the resulting response function is then confronted to 
what one obtains using the ESRPA by performing a perturbative expansion of the responses in 
each case. The possibility of having a three nucleon ejection process is finally analyzed in this 
framework. 

Let us start with the linear response to an external field F defined as 



R{E) = -I / (0|r F^t(^)^^(o) |0)e'^*dt. 



(1) 



where F (t) = e Fe , H = Hq + V, with Hq and V being, respectively, the Hartree- 
Fock (HF) mean field and the residual interaction. The spectral representation of the response 



function, in terms of a set {\iy)} of eigenstates of the hamiltonian H, reads 

■(0|F|z/)(z/|Ft|0) (0|Ft|z/)(z/|F|0)' 



RiE) = Y. 



E-E^ + ir] E + E^-ir] 

where rj is an infinitesimal positive number. 

Within the equation of motion method m, the set {|z^)} is generated as 



J ^v 



and 



(2) 



(3) 



(4) 



n^\0) = 0, for all u. 

The operators CJ and Ci (with CJ = aj,^ ■ ■ ■ al,.ah^ ■ ■ ■ ah^) create and annihilate i particle-hole 
pairs on the HF vacuum |0) = \OpOh), respectively. 
The equation of motion for Q^ 



a. 



H,ni 






|6)5.,M, 



(5) 



|0) = E.(0| 

where E,y stands for the excitation energy of the state \i>), leads to the RPA-like eigenvalue 
problem 

AX'' = EMX", (6) 

with 



The submatrices A, B and A^ given by 






M 



N 

-A^* 



(7) 



A 



«j 



a 



H, C]]] |0), B,, = (0| [Q, [H, C,]] |0), N,, = (0| [Q, C}\ |0), (8) 

and using eqs. (H-|D it is possible to write the response in representation independent form as 

R{E) = T\EJ\r -A + iriiy^T, (9) 

where T is defined as 



^ =\ pB ) ' with 



{o\[a,F\\o), 

pB ^ pM^p _^ p]y 



(10) 



After splitting the Hilbert space of ipih states into a P-space that includes only the Iplh 
states and the Q-space that spans on the rest of the states, the response function can be written 
as 



where 



with 



and 



R{E) = Pp{E)Gp{E)Tp{E) + mQ{E)TQ, (11) 



gp{E) = [EAfp + ir]Ip -Ap- {ApQ - AfpgE) QgiE) {Aqp - MqpE)]-' , (12) 



QgiE) = [EMq + i^Xq - AqY^ , (13) 



^p{E) =Tp- MpqTq + ApqQq{E)Tq. (14) 



In standard RPA the state |0) is approximated by the HF ground state and the Q-space 
is absent, while the so called extended RPA incorporates perturbative ground state 2p2/i ad- 
mixtures and a perturbatively suggested truncation of the dynamical matrices and excitation 
operator. It is obtained by: 
i) evaluating the matrix elements (|^) and ( p!0D for |jll| 



|0)=co|0) + ^C2o|2o), (15) 

2o 

where 

C0 = 1-^E|C20I% C2o = -^, (16) 

^ 2o ^2o 

2o = (piP2^i^2)o represents the 2p2/i ground state admixtures, i?2o the corresponding unper- 
turbed energy and V200 = (2o|^|0), and 

ii) keeping terms up to second order in V for the forward sector within the P space, terms 
linear in V for the backward sector within the P space and for the coupling between the P and 
Q spaces, and only terms of zeroth order within the Q space. Under these conditions the norm 
matrix elements read 0) 

Ni^ = 5,j + AN, J (17) 



where i = ipih and the nonzero ANij are 



ANw = E 4,C2',{2o\Dii'\2'o), AN,s = ^^^^(l; 2o|3), (18) 

20,2;, 2o 

where Dw = Ci, C\i — dw and (1; 2o|3) is the overlap between the Iplh ® {2p2h)o and 3p3h 
final state configurations. (Note that within the quasi-boson approximation Dur = 0.) The 
explicit result for the matrix element (2o|-Dii'|2q) is 



{{piP2hh)o\D,i,^p'h'\{p[P2KK)o) = - [1 + P{hu h2)P{h[, h'^)] 
Sp,p'Shi^h'P~{h, h2)P~ {pi, P2)Sh[,hSh2,h'2^P2,p'2^Pi,p[ +P "^ h, 



(19) 



where P {i,j) = [1 — P{i,j)], while the operator P{i,j) exchanges the arguments i and j. 
The forward going energy matrix elements are evaluated in the same way and one gets 



Aij = 5ijEj + Vij + AAij, 
where Vij = (^|V^|j) and the nonzero matrix elements AAij are: 

AAn'= E(^i-^2o)c;„C2^(2o|Ai'|2[,), AA,s = AN.^E^. 

2o,2[, 

The one-body matrix elements are: 

fi + Ei'ANu'fi' hTt = l 

E20 C2o/i2o for i > 1, 

where 

/i = (l|F|0), and /,2„ = (^|F|2o). 

Before proceeding it is convenient to introduce the unperturbed Green's function: 



(20) 



(21) 



F/- 



Q\E) 



G\E) 

G'°*(-E) 



(22) 



(23) 



(24) 



where G^{E) = [E^ — A{H = Hq)] ^ (with E^ ^ E + irf) and rewrite the perturbed Green 
function within the space P in the form 

Qp{E) 



where 



with 



and 



lCp{E) = lC^v{E) 



{gliE))-' - lCp{E) 



Bill 






(25) 
(26) 
(27) 



i=2,3 



^(2) 



AJ:\\i{E) = AAir-ANirE, 
Asg(E) = 



2Vis-ANi,(Gl,{E)) )AN,r. 



-1 



(2^ 



In the above equations V,y stands for the matrix representation of the residual interaction within 
the ipih ® jpjh subspace. 

The response function now reads 



R{E) = Tv{E)gii,{E)h,{E) + Yl J'!Gu{E)y'u 

i=2,3 



(29) 



where 



:FiiE) ^ 



F^E) 
Fi^iE) 



, with < 



^ F,^{E)=fi + AFiiE), 
AFi{E) = AFf + AFf + Ei=2,3 Vi;G'^XE)F; (30) 

I AFf = AiVn'/i', AFf = -AN^Fs. 

From the expressions for AAw and ANu/, given by Eqs. (|T8|) and ( pi]) respectively, the 
matrix elements AS![]^, (F) and AF| can be expressed as: 

AeS(^) = -j:^*2o^%i^o\Dn'\2'o)iE-Ei + E,,), (31) 

2o,2^ 

AFf) = ^ c;„C2;, (2o|Ai'|2[,)/i'. (32) 

2o,2(, 



Moreover, from the relationships 



^i3 = -Ec;o^2o(l;2o|3); /3,2o=E(3|l;2o)/i, (33) 

2o 1 

one obtains 

AEg = ^c;„C2^(l;2o|l';2[,)(E-Ei + E2i), (34) 

2o,2^ 

AFf) = -^c;„C2^(l;2o|l';2;,)A,. (35) 

2o,2(, 

We can note here that 

(1; 2o|l'; 2[,) = (2o|(Ai' + div)\2'^), with Jn, = Sir + ^Ci, (36) 

and thus in summary we get: 

i) in the ESRPA (where the Q space includes only the 2p2h excitations) 

Sn'(i?) = -j:^*2o^2',{2o\Du'\2',)iE-Ei + E,>^)+J2^^^r, (37) 

2o,2i 2 E+-E2 

UE) = h+ Y: c;„C2;,(2o|^ir|2-)/r + E ^'/!"!?" ; (38) 

2o,2;,;l' 2,2o; ^ "^2 

ii) in the ETRPA (where the Q space includes both the 2p2h and 3p3h excitations) 

^n'iE) = Y^ clc,,{2o\du'\2',)iE - E, + E,,) + ^ ^^, (39) 

20,2;, i=2,3 -^ -^* 

2o,2j, i=2,3;2o ^ ^^ 



The results (|37|) and (|38D are in essence those obtained previously by Arima and collaborators 



P, |TT| and by the Jiilich group p, |T^. On the other hand, when terms containing the matrix 

elements {2o|(iii'|2Q) are neglected in eqs. ( p9D and (^0[), one finds the results derived in our 

previous works [0]. Q 

^These terms give rise to disconnected graphs, which are nonphysical, as weU as to double connected graphs 
represented in fig. Id and le, respectively. As seen from relations (E3) and ( p7| ) below, they do not contribute 
to the response function. 



In order to elucidate some of the content of these equations we turn next to a perturbative 
expansion of the response function and examine the leading corrections to the unperturbed 
Iplh response R^{E) = X^i |/iP/(-^^ " -^i)- To achieve maximum simplicity we first omit the 
residual interaction within the Iplh sector and backward contributions, so that to second order 
the Bethe-Salpeter equation Eq. (|25|) reads 



Gu'iE) = Gl,{E) + Gl,{E)J:n'{E)G'i,,,{E), 



(41) 



which substituted in Eq. (pH]) leads to the desired approximation for the response function. 
Within the ESRPA one gets: 

R{E)=R{E)+ ^ c,„^;^— ^C2^+2^>_^^-^^— ^^^— ^ 



2,2o,2[, 



'E+-Eo 



+ E' 



/r 



1,1' ^~^ ~ ^^ 
and in the ETRPA: 



2o,2{, 2 



E+-E,. 



h 



E+ - Er 



(42) 



i?(E) = /?°(E) + E C2o#%c,,+2 E ^^^^*^^°'^°^ ^'' 



J=2,3;2o,2^ 



'E+-E, 



J=2,3;l,2o 



E+ -El E+ - Ei 



n 



y 



Y clc,>^{2o\duWiE -El- E,,) - Y 
20,2;, 



1=2,3 



^^13^31' 

E+ - E, 



h 



E+ - Ey 



(43) 



Now the two expressions (|42|) and (^31) can be shown to be equivalent. This results in fact 
from explicitly performing the sums over 3p3/z. states in Eq. (143|) . To do that one first rewrites 
these sums making use of relations (^3]) and (|36|) as: 



E 0. 

3,2o,2^ 



* /3*2o/32i ^ Y^ ^* ^* (2o|(-Pll' + C^ll')|2o) f ^ 
^2(, = Z^ C2J1 r;:^ /I'Ca;, 



■^0 E+ - E-i 



3,2o,2(, 



£,+ — ]^Y — _£,, 



(44) 



2o 



2E 



3?(/l732oC2o) ^1 



13 



l,3,2o 



E+ -El E+ - E3 



. >r^ * ,. i?2o(2o|(Ai^+rf;iO|2;,) ^ 
"2 2^ ClJlTT^ FTTT^T ^ Fr^/l'C2;, 



l,2o,2[, 



{E+ - Ei){E+ - El - E. 



2oy 



(45) 



and 



^^^E+-E^E+-E^E+-Ev ^^^f^^^^,^{E+-Ei){E+-Ei-E2,){E+-Eyy ^ ^ 



The result of performing the sum is: 

1,1 ,2o,-iQ ^ 



(47) 



which substituted in Eq. (|43| ) gives the expression ( ^21) also for the ETRPA response. 

The cancellation among the 3p3/i on the energy-shell contributions can be exhibited also 
making use of the Rayleight-Schrodinger perturbation expansion, i.e., 

\i) = \i) + \i)W + . . . and Ei = Ei + ^f ^ + ■ ■ -, i = ipih, (48) 

where the perturbed wave functions and energies are indicated by the symbol ~ and the su- 
perscript points the order of the correction introduced by the residual interaction V on the 
unperturbed quantities \i) and Ei. The amplitude for the E- excitation from the correlated 
ground state to the perturbed 3p3h states reads 

mo)^^Mm.m«JM^oiv% (49) 

£/3 - Eq rys - J^o 

with 

(3|[^, F]|0) = E(3|V^|1)(1|^|0) - E(3|^|2o)(2o|V^|0) = 0, (50) 

1 2o 

where the last equivalence is a direct consequence of the relations (p3D, i.e., 

E^3l/l = -Ec2oi?2o(l;2o|3)/l = E/3,2oV2oO (51) 

1 l,2o 2o 

Thus we see once more that, up to the second order in V, the 3p3h final states do not 
contribute to the response function and that |(3|F|0)p = OiV"^). The Goldstone diagrams for 



^Note that Ho does not contribute since {3\[Ho,F]\0) ~ 0. 



the fourth order 3p3h on the mass-shell contributions to the response function are shown in 
figs, le and If. 

At first glance it might look as if the connected Goldstone diagrams associated with the terms 
(@1)' (HD ^^^ (HI) of ^he ETRPA response (illustrated in Figs, la, lb and Ic, respectively) 
should give rise to on the mass-shell 3p3h contributions, through the imaginary part of the 
propagator {E^ — E's)"^. However, Eq. (^) shows that these contributions in fact cancel out 
so that the 3p3h sector only affects the Iplh excitations by coupling them with the virtual 
intermediate states |l;2o). Thus in spite of including the 3p3h propagator in the Green's 
function, three nucleon ejection does not occur in the leading order processes. The above 
mentioned diagrams also explain the physical meaning of the fourth term in the expression 
(|42|). The cancellation of on shell 3p3h contributions results from the destructive interference 
between amplitudes involving creation of the 3p3h state from a ground state correlation and 
from V31 coupling respectively. A similar calculation in which the backward part of Eq. ( pSf ) 
and/or the residual interaction within the Iplh space are kept up to the relevant order leads 
again to the same result. It is worth stressing that this does not depend on the form of the 
two-body force used as residual interaction or on the size of single particle space. 
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Figure 1: Graphical representation of the second and fourth order contributions to the response 
function. The dotted circles (©) denote the one-body vertices and the filled ones (•) indicate 
the two-body matrix elements. The diagrams (a), (b) and (c) correspond, respectively, to the 
terms given by eqs. (il), (|45| ) and (^Hf ). Second order unlinked and double-linked graphs 
analogous to the diagram (c) are shown in figures (d) and (e), respectively. The last ones, 
although contained in eqs. (|39D and (^OD , do not contribute to the response function. Finally, 
figure (f) illustrates the fourth order on the energy-shell 3p3h processes. 
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